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We present electronic transport and magnetic properties of single crystals of 
semimetallic half-Heusler phase LuPdBi, having theoretically predicted band 
inversion requisite for nontrivial topological properties. The compound exhibits 
superconductivity below a critical temperature = 1.8 K, with a zero-temperature 
upper critical held Bc 2 ~ 2.3 T. Although superconducting state is clearly rehected 
in the electrical resistivity and magnetic susceptibility data, no corresponding 
anomaly can be seen in the specihc heat. Temperature dependence of the electrical 
resistivity suggests existence of two parallel conduction channels: metallic and 
semiconducting, with the latter making negligible contribution at low temperatures. 
The magnetoresistance is huge and clearly shows a weak antilocalization effect in 
small magnetic helds. Above about 1.5 T, the magnetoresistance becomes linear 
and does not saturate in helds up to 9 T. The linear magnetoresistance is observed 
up to room temperature. Below 10 K, it is accompanied by Shubnikov-de Haas 
oscillations. Their analysis reveals charge carriers with ehective mass of 0.06 rrie 
and a Berry phase very close to tt, expected for Dirac-fermion surface states, thus 
corroborating topological nature of the material. 
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Large family of half-Heusler compounds that crystallize in a noncentrosymmetric cu¬ 
bic MgAgAs-type structure attracts much attention due to their remarkable magnetic and 
electrical transport properties.^ Recent ab initio electronic structure calculations have pre¬ 
dicted that several dozen of half-Heusler compounds, due to a strong spin-orbit coupling, 
have inverted bands requisite for topological properties.®^ Topological insulators (TI) con¬ 
stitute a new class of materials, which are insulating in the bulk, but at the same time their 
surface states are protected from backscattering by Z 2 topology.®® These surface states are 
theoretically described as massless Dirac fermions with linear dispersion and lifted spin- 
degeneracy. Such properties open new prospects for spintronic applications. In contrast to 
two-dimensional (2D) TIs, like HgTe, some of half-Heusler phases possess not only band 
inversion at an odd number of time-reversal-invariant momenta, but also may have a bulk 
band-gap opening under uniaxial pressure, which lowers the crystal symmetry.® In result, 
the rare-earth {RE) based half-Heusler compounds (RARdBi, RARtBi, Ri^PdSb, Ri^PtSb) 
may possibly form the biggest group of the three-dimensional (3D) TIs. 

After the discovery of the 3D TI systems, a quest for superconductivity in materials 
with the non-trivial topology of electron bands started due to their potential applications in 
topological quantum computing.^I®!^ Topological superconductor (TS) is a material charac¬ 
terized by protected Majorana surface states and bulk consisting of mixed-parity Cooper-pair 
state^i^ni! The hrst theoretically forecasted and experimentally discovered TS was Bi 2 Se 3 
doped with Cu atomd^^^. Since it remains problematic to clarify the nature of the surface 
Majorana fermions in the inhomogeneous crystals of Cua:Bi 2 Se 3 fPef. ITTj) . new TS materials 
are sought extensively among low-carrier-density semiconductors, whose Fermi surfaces are 
centered around time-reversal-invariant momenta.^ In this context, the superconductivity 
recently found in the bismuthides REPtBi: LaPtBi {T^ = 0.9 K),^ LuPtBi (Tc = 1.0 K),®' 
and YPtBi {R = 0.77appears most attracting. Besides, the content of magnetic 
rare-earth may bring about other properties, such as long-range antiferromagnetic ordering 
or heavy fermion behavior.l^ 

Potentially, the rare-earth palladium-bearing bismuthides Ri^PdBi are equally interest¬ 
ing from the point of view of topological and superconducting properties. Magnetic and 
transport behaviors of these compounds have been hrst investigated on polycrystalline 
samples.l^^®®' Most of them have been characterized as local-moment antiferromagnets 
(with Neel temperatures Tn = 2 — 13 K) with electrical properties characteristic of semimet- 
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als or narrow-band semiconductors. Recently, diamagnetic YPdBi has been reported to 
exhibit large linear magnetoresistance (LMR) and Shubnikov-de Haas (SdH) quantum oscil¬ 
lations at low temperatures, and its nontrivial topological nature has been evaluated.^ In 
turn, CePdBi has shown some features of incipient superconductivity (Tc = 1.3 K) emerging 
in the antiferromagnetic state (Tn = 2K), and a hypothesis on topological character of the 
compound has also been formulated.!^ Similarly, ErPdBi has been reported to show the 
coexistence of antiferromagnetism (Tn = 1.06 K) and superconductivity (Tc = 1.22 K), and 
hence acclaimed as a new platform for the research of the reciprocity of magnetic order, 
superconductivity, and nontrivial topological states.!^ However, in our own study on single¬ 
crystalline ErPdBi, we could not conhrm those Endings because no clear-cut evidence for 
any intrinsic superconducting state in this material has been observed.!^ Most recently, su¬ 
perconductivity (Tc = 1.7 K) has been reported for nonmagnetic LuPdBi, accompanied by a 
2D weak antilocalization (WAL) efiect.!^ Since the latter feature is considered as a fingerprint 
of topological surface states, the compound was suggested to be a topological superconduc¬ 
tor with Majorana edge states.!^ In this work, we present the results of our comprehensive 
magnetic susceptibility, specific heat and electrical resistivity measurements performed on 
high-quality single crystals of LuPdBi. The obtained data are here critically compared with 
those reported in Ref. [3l 


RESULTS 

Electrical resistivity and Hall effect 

The temperature dependence of the electrical resistivity, p, of the LuPdBi single crystal is 
shown in Figure [T^. From 300 K down to ~ 170 K, the resistivity displays a semiconducting- 
like character (dp(T)/dT < 0), then at lower temperatures p{T) becomes typical for metals. 
Such a behavior of the resistivity is common for the RTPdBi compounds.®®^^^ As demon¬ 
strated in Figure [^, the overall behavior of the electrical resistivity in LuPdBi can be well 
described over an extended range of temperature considering two independent channels for 
charge transport: semiconducting- and metallic-like. The former contribution is written as 
PsiT) = aexp(—T^/Z cbU), where Eg is an energy gap between valence and conduction bands, 
while the latter one is expressed as PmiT) = po + bT"^ + cT, where po is the residual resis- 
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FIG. 1. Electrical transport properties of single-crystalline LuPdBi. (a) Temperature 
dependence of the electrical resistivity up to 300 K. Inset: low-temperature resistivity data revealing 
a superconducting transition near 2 K. (b) Temperature dependence of the electrical conductivity in 
the normal state fitted by the function cr{T) = ara{T) + as{T) (red line) described in the text. The 
metallic and semiconducting contributions are represented by blue and green curves, respectively. 

tivity due to scattering on structural defects, the second term represents electron-electron 
scattering processes, while the third one accounts for scattering on phonons. The total con¬ 
ductance, o'(T) = l/p(T) is the sum cr(T) = as(T) -|- crni(T), where o's(T) = l/ps{T) and 
(Tm(T) = l/pm(T). Fitting these expressions to the experimental data of LuPdBi in the 
wide temperature interval from 5 K to 300 K, yielded the parameters: a = 296 mGcm, Eg 
= 11.5 meV, Pq = 226mf2cm, b = 0.018 mf2cm/K^, and c = 1.72mGcm/K. The so-derived 
energy gap Eg is very close to 15.2 meV reported for LuPdBi in Ref. [TUI very similar to 
the values reported for other RF'PdBi phases.®®^!*^ 

It is worth noting that the temperature variation of the resistivity of LuPdBi single crystal 
obtained in the present work somewhat differs from that reported before by Xu et al.^, who 
observed a monotonous increase of the resistivity with decreasing temperature from 300 to 
2 K, and attempted to analyze their p{T) data above 90 K in terms of 3D variable-range 
hoping model. However, also in that case, the transport at low-temperatures was attributed 
to the metallic channel (though no quantitative analysis was made). 

As displayed in Figure [^, the Hall resistivity of LuPdBi measured at 2.5 K is a linear 
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FIG. 2. Hall effect in LuPdBi. (a) Magnetic field dependence of the Hall resistivity measured 
at 2.5 and 300 K. (b) Temperature variation of the Hall carrier concentration, (c) Temperature 
variation of the Hall carrier mobility. 


function of the applied magnetic field up to 9 T. The Hall coefficient evaluated from this 
data was Ru = 0.0565 mGcm/T. Almost identical magnitude of Ru was derived from Pxy(-B) 
collected at T = 4, 7, 10 and 50 K (not shown). Assuming a single parabolic band, one may 
estimate the Hall carrier concentration, uh, to be at low temperatures uh — 1-2 x lO^^crn”^. 
Above 50 K, uh was found to increase with increasing temperature and to attain at 300 K 
a value of 2.8 x 10^®cm“^ (see Figure]^). In turn, as demonstrated in Figure]^, the so- 
obtained Hall carrier mobility, pn, gradually decreases with increasing temperature, from 
the value 2404cm^V“^s“^ at T = 2.5 K to 573cm^V“^s“^ at 300 K. 

The inset to Figure [pi displays the resistivity of single-crystalline LuPdBi at the lowest 
temperature studied. In agreement with Ref. |3l the compound undergoes a transition to 
superconducting state. The resistivity starts to drop near 2.2 K and becomes zero below 
1.8 K. Using a criterion of 90% drop in the resistivity value from the point where the 
p{T) curve starts bending down, T^. = 1.9 K was determined. This value is slightly higher 
than Tc = 1.7 K reported before for LuPdBi (Ref. |3j) and remains the highest R found 
so far for the RETBi superconductors. As can be inferred from the inset to Figure 
besides the superconducting phase transition, we observed a small resistivity drop below 
5 K. The existence of a fairly similar anomaly in p{T) near this temperature has previously 
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FIG. 3. Superconductivity in LuPdBi. (a) Low-temperature dependencies of the electrical 
resistivity measured in different external magnetic fields, (b) Upper critical field as a function of 
temperature. Solid line represents the fit of the function Bc 2 {T) = B*2{1 — T/Tc)^~^‘^ described in 
the text. 

been reported for a few other half-Heusler bismuthides and antimonides,®® and tentatively 
attributed to superconductivity of thin hlms of metallic Bi(Sb) located at grain boundaries, 
yet the actual nature of this feature remains unclear. 

The temperature dependence of the electrical resistivity of LuPdBi measured in external 
magnetic helds up to 2 T is shown in Figure [^. The width of superconducting transition 
increases gradually from 0.4 K at zero held to 0.7 K at 1 T. The critical temperature dehned 
by the drop of the resistivity to 90% of its normal-state value varies with magnetic held in a 
manner presented in Figure]^. The so-derived Bc 2 {T) dependence has a positive curvature, 
characteristic of two-band clean-limit type-II superconductors, like YNi 2 B 2 C, LuNi 2 B 2 C 
(Ref. ESI) MgB 2 (Ref. ES])- Though our sample exhibits the residual resistivity ratio 
p(300 K) /p(4 K) of only 1.8, the observation of SdH oscillations described below, indicates its 
very good quality, which allows us to consider the observed superconductivity to have a clean 
limit character. Indeed, as displayed in Figure]^, the experimental Bc 2 {T) data can be well 
approximated in the whole temperature range by the expression Bc 2 {T) = B* 2{1 — T, 
previously applied, e.g., to MgB 2 (Ref. ESI- The htting parameter B *2 = 2.3 T can be 
considered as the upper limit for the upper critical held i?c2(0). The parameter a = 0.496(22) 
is signihcantly larger than 0.32 and 0.24 (0.25) reported for MgB 2 and LuNi 2 B 2 C (YNi 2 B 2 C), 
respectively.^^^® Using the relation ,^o = ^*ho/2vri?c2(0), one obtains an estimate for the 
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coherence length = 17 nni. This value is about three times larger than those reported 
for MgB 2 (Ref. EHl), YNi 2 B 2 C and LuNi 2 B 2 C (Ref. EZj). Nevertheless, it is still distinctly 
smaller than the mean free path I = 60 nm derived from the SdH data (see below) and this 
comparison corroborates our assumption that the superconductivity in LuPdBi occurs in 
the clean limit {l/^o > !)• The Werthamer-Helfand-Hohenberg approximation!^ provides an 
estimate for the orbitally limited upper critical held Rorb = 0.72Tc[—dRc 2 /dT]T^ to be 0.5 T. 
The Pauli limited upper critical held Bp can be evaluated from the equation Bp = A/\/2ftp, 
where the BCS energy gap A = 1.76 kpTc. The estimate for the measured crystal of LuPdBi 
is Bp =3.5 T. Thus, the Maki parameter ctm = Borh/Bp = 0.2 was obtained. One can 
conclude that the superconductivity in LuPdBi is Pauli limited. It is worthwhile recalling 
that the relationship Rorb < Rc 2 < Rp was established also for the related half-Heusler 
superconductors YPtBi and LuPtBi (Refs. EO] and |22)) . 


Thermodynamic properties 

Figure [^displays the results of magnetic measurements performed on a collection of sev¬ 
eral single crystals of LuPdBi. Both components of the AC magnetic susceptibility show 
clear anomalies at 1.8 K typical for a superconducting transition. The corresponding feature 
is seen on the temperature dependence of the DC magnetic susceptibility measured in a weak 
magnetic held upon cooling the sample in zero held. Comparison with the data obtained 
in held-cooled regime reveals superconducting current screening ehect. Furthermore, the 
magnetization M measured at the lowest temperature T =1.62 K attainable in our exper¬ 
imental setup exhibits a held variation characteristic of type-II superconductors. All these 
features are in concert with the electrical resistivity data and manifest the emergence of the 
superconducting state in LuPdBi below Tc = 1.8(1) K. From the M{B) curve presented in 
Figure one can roughly estimate the upper limit for the lower critical held R^ to be of 
0.05 mT. Then, the formula Rd = d>o/27rAg(, yields the penetration depth Asc = 3629 nm. 

As can be inferred from Figure the low-temperature dependence of the specihc heat 
of LuPdBi is featureless. In particular, no obvious anomaly is seen at the superconducting 
phase transition. Remarkably, analogous situation occurs for other half-Heusler supercon¬ 
ductors: LuPtBi (Ref. [Tj) and YPtBi (Ref. |8]). Lack of any feature in C'(T) at Tc in these 
systems may be due to the fact that superconductivity originates from topologically pro- 
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FIG. 4. Magnetic properties of single-crystalline LuPdBi. (a) Temperature variations of 
the real and imaginary components of the AC magnetic susceptibility taken with a magnetic field of 
2.7 /xT alternating with a frequency of 317 Hz. (b) Field dependence of the magnetization measured 
at T = 1.62 K upon cooling the specimen in zero magnetic field, (c) Low-temperature dependence 
of the DC magnetic susceptibility recorded upon cooling the specimen in zero magnetic field (ZFC) 
and in an applied field of 40 /xT (FC). All solid lines are guides to the eye. 

tected surface states, which occupy a very small volume in relation to the volume of entire 
sample. The only exception is the result reported by Xu et al.^who observed for their single 
crystal of LuPdBi a prominent BCS-like specific heat jump at the onset of superconductivity. 
A critical comment on this data can be found in Supplementary Material. 

As shown in the inset to Figure]^ below 5 K, the C/T{T) variation can be well described 
by the standard formula C/T = 7 -t- with the electronic specihc heat coefficient 

7 = 0.75mJ/molK^ and the coefficients of the phonon contribution jS = O.SmJ/molK"^ and 
5 = 0.02 mJ/mol K^. From the relation ©d = ( 12 ni 77 r^/ 5 / 3 )^/^ {R is the gas constant and n is 
the number of atoms per formula unit) one derives the Debye temperature ©d = 194 K. The 
obtained values of 7 and ©d are close to those reported in Ref. ITOl 7 = 0.096(327) mJ/mol 
and ©D = 194.9 K but strongly disagree with the results derived for LuPdBi by Xu et 
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FIG. 5. Heat capacity of LuPdBi. Low-temperature dependence of the specific heat. Inset: 
specihc heat over temperature ratio as a function of temperature squared. Solid red line represents 
the ht to the function CjT = 7 -|- /ST^ -|- discussed in the text. 

al.: 7 = 11.9mJ/molK^ and ©d = 295K (Ref. [3]). It is worth recalling that for other 
nonmagnetic RETJii the valnes of 7 are also very small: 0.1 for YPtBi (Ref. HT|) and 
0.06mJ/molK^ LnPtBi (Ref. [7]). 


Magnetoresistance 

Fignre shows the magnetoresistance data of LnPdBi, defined as MR(R) = [p{B) — 
p(0)]/p(0)], with p{B) and p(0) denoting the electrical resistivity measnred in applied field 
and in zero held, respectively. Remarkably, at 2.5 K MR is as large as 210% in the strongest 
held of 9 T achievable in onr experimental setnp. With increasing temperatnre, MR gradn- 
ally diminishes, yet even at 300 K it reaches abont 30% in R = 9 T, which is a very large 
valne. A striking featnre of MR(R) is its linear nonsatnrating behavior. Above a certain 
magnetic held (ranging from 1.5 T at 2.5 K to 3 T at 300 K), MR becomes proportional to 
B with a slope systematically decreasing with increasing temperatnre, from abont 8 %/T at 
the lowest temperatnres, down to 3%/T at 300 K. 

In conventional theory of metals and semicondnctors, MR shonld grow qnadratically in 
weak magnetic helds and then satnrate in strong helds.®' Linear magnetoresistance (LMR) 
is a phenomenon reported for narrow-gap semicondnctors,l22H15l mnlti-layered graphene,®! 
topological insnlatorP^HlSI ^^d Dirac semimetals.®!^ Two main approaches were made to 
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FIG. 6. Transverse magnetotransport in single-crystalline LuPdBi. (a) Field dependence 
of the magnetoresistance measured at different temperatures. Solid lines emphasize the linear 
behavior, (b) Low-field variations of the conductance measured at different temperatures. Solid 
curves represent the functions given in Eq. which account for the weak antilocalization effect. 

understand this behavior. In their classical theory, Parish and Littlewood attributed LMR 
to material inhomogeneities, which give rise to sturdy spatial fluctuations in the electrical 
conductivity.!®' In turn, Abrikosov developed the quantum magnetoresistance theory,!®! that 
predicts LMR in zero-gap band systems with linear energy dispersion being in the extreme 
quantum limit (EQL) in which all electrons are confined to the first Landau level. The 
Abrikosov theory was demonstrated to work correctly for graphend®! and topological insula¬ 
tor Bi 2 Se 3 nanosheets,!®! with topological surface states possessing a linear energy dispersion 
over wide temperature interval. 

For single-crystalline LuPdBi, the detailed X-ray and microprobe characterization re¬ 
vealed a very good sample quality, which was further confirmed by the observation of SdH 
oscillations (see below). This prompts us to dismiss the classical theory. On the other hand, 
the compound shows LMR behavior virtually identical to that reported for silver chalco- 
genides, interpreted within the quantum magnetoresistance theory.!®! Moreover, the SdH 
data indicated that the studied crystal was in an EQL regime at temperatures up to at least 
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10 K. All these findings suggest that the quantum approach is perfectly applicable, and the 
observed LMR effect is a hallmark of Dirac-like band structure in this putative topological 
semimetah 

In weak magnetic fields, MR of LuPdBi displays a very abrupt rise. For example, at 
T = 2.5 K, MR attains a huge value of 150% already in a field of 1.5 T. Similar situation, 
but with lower MR increase, occurs at higher temperatures, up to 150 K. Such sharp cusps 
on the MR isotherms are reminiscent of a weak antilocalization (WAL) effect, which relies on 
destructive interference of the electron wave functions between two closed paths with time- 
reversal symmetry. In topologically nontrivial materials, WAL is enhanced due to ability 
of topologically protected surface states to collect a vr Berry phase due to the helical spin 
polarization. Thus, the WAL effect may be considered as a manifestation of surface-state 
transport in LuPdBi. 

Taking into account the coexistence in our LuPdBi sample of two parallel conduct¬ 
ing channels (see Figure and the discussion above), in conjunction with the plausibly 
topological nature of the compound, one may assume that at temperatures below ^ 25 K 
the conductivity was almost entirely determined by the metallic channel due to the sam¬ 
ple surface. Consequently, the sheet resistance Rs{B) of the conducting surface of the 
sample can be derived. Figure displays the sheet magnetoconductance of LuPdBi, 
Aa{B) = {l/Rs{B)) — (l/i?s(0)), where Rs{B) and Rs{0) represent the data obtained 
in applied field and zero field, respectively. These experimental results can be analyzed in 
terms of the Hikami-Larkin-Nagaoka formulaP^ 

- R h \ , r h 


Aa{B) = 


rje 
2n‘^h - 


2 ^ AeLlB 


ln(^ 


AeLlBR^ 


( 1 ) 


where stands for a phase coherence length, r; is a parameter depending on the strength 
of spin-orbit interaction and magnetic scattering, 'ip{x) is digamma function, while the 
other symbols have their usual meanings. For T = 2.5 K, the fitted parameters are 
Lp = 49.6(2.2) nm and rje^/{2'K‘^fi) = —153(21) and are comparable to Lp = 95 nm and 
rje^/{2Ti‘^h) ~ —10 obtained for the same compound in Ref. O For higher temperatures, 
where the contribution due to the bulk channel becomes noticeable, the fits of Eq. [T] to 
the experimental data are still of good quality, but the value of r] decreases with increasing 
temperature. More detailed description of the dependence of WAL parameters on temper¬ 
ature is given in Supplementary Material. It is worth mentioning that a similar analysis of 
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FIG. 7. Shubnikov—de Haas oscillations in LuPdBi. (a) Oscillatory component of the elec¬ 
trical resistivity as a function of inverse magnetic field, measured at different temperatures. Black 
solid line represents a fit of 2.5 K data with L-K formula. Inset: FFT spectra obtained for different 
temperatures, (b) Oscillatory component of longitudinal electrical conductivity at T = 2.5 K cal¬ 
culated with macroscopic Hall resistivity (blue line) or surface channel Hall resistivity (red line). 


MK{B) of LuPdBi was presented in Ref. |3l however the fitted values of r] were ~ 10^ larger 
than that expected for a TI system (see also our comment in Supplementary Material). 


Shubnikov-de Haas oscillations 

SdH oscillations are being recognized as a powerful tool in the studies of 3D TIs.^ Dirac 
fermions are characterized by a tt Berry phase (for ordinary metals the Berry phase is equal 
to zero), and the phase factors of the oscillations reveal the Berry phase values for the 
system. Observations of SdH oscillations have been reported in some candidate topological 
materials from the half-Heusler family. ^P^P^P^ * For our LuPdBi samples, SdH oscillations 
were well resolved at temperatures up to 10 K, hence indicating very good quality of the 
single crystals measured. Figure shows the oscillatory component Apxx, obtained by 
subtracting a background from the Pxx{B) isotherms. The oscillations are clearly periodic 
in 1/R, however their frequency is very low. Fast Fourier transform (FFT) analysis of 
these data yielded rather poor quality spectra with broad peaks corresponding to a single 
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temperature-independent frequency ~ 10 T, shown in inset to Figure igi. Obviously, the 
broadness of FTT peaks comes from low frequency of the oscillations and only less than 
three periods covered by our measurements. In order to obtain more reliable value of the 
oscillation frequency /sdu we perfomed direct htting of the standard Lifshitz-Kosevich (L-K) 
expression where Apxx ~ cos[27r('^^ — V 2 + /d)] with (3 being a phase shift The analysis 
yielded /sdu = 11-9 T and {3 = 0.431(2). 

The oscillation frequency is related to the extremal cross section area A-p of the Fermi 
surface through the Onsager relation /sdu = {h/2'Ke)Ap, where Ap = nkp. The value of 
/sdH derived for LuPdBi implies the Fermi vector kp = 0.019 that corresponds to a 2D 
carrier density ^20 = kp/An = 2.87 x 10^^ cm“^ (we assume lifted spin degeneracy). The 
temperature dependence of the resistivity oscillations amplitude is given by the standard 
L-K expression Apx^{T) ~ A(T)/sinh(A(T)), where A(T) = 27r‘^k-BTm*/heB and m* is the 
effective cyclotron mass.^ Fitting this equation to the experimental data of LuPdBi (see 
Figure]^) one obtains m* = 0.06me, where me is the free electron mass. From the values 
of kp and m*, the Fermi velocity up = hkp/m* ~ 3.6 x 10^ m/s and the Fermi energy Ep = 
m*Vp ~ 45.7 meV were derived. Since calculation of the Sommerfeld coefficient for a 2D 
Fermi surface is impossible (it encloses zero volume), we used the formula 7 = k‘^Vm*kp/?th^ 
(corresponding to spherical Fermi surface) and obtained the Sommerfeld coefficient equal to 
0.003 mJ/mol K^. This value is much smaller than 0.75mJ/molK^ derived from the specihc 
heat data. This hnding indicates that the SdH oscillations originate from 2D states Fermi 
surface or from a tiny 3D Fermi surface, much smaller than the bulk one contributing to the 
specihc heat. If the bulk states were contributing to SdH oscillations one should observe high 
cyclotron frequency (large kp) and/or large cyclotron mass, not seen in LuPdBi. Moreover, 
the Hall concentration of uh — 1.2 x 10^®cm“^ is over 50 times larger than 2.3 x lO^^cm”^ 
obtained when we assumed that SdH effect originates from 3D spherical Fermi surface. Such 
discrepancy in the Hall and SdH concentrations further supports the conjecture that the SdH 
oscillations in LuPdBi originate from the surface states. 

Knowing m*, one can carry out the Dingle analysis shown in the inset to Figure [^. The 
lifetime r of the carriers can be found from the Dingle temperature Td = h/2nkBT. The slope 
of the linear £t to the data yielded Td = 16.9 K, which gave r = 7.2 x 10“^^ s. Ergo, the mean- 
free path I = vpT = 26 nm and the surface mobility ps = er/m* = el/hkp = 2100 cm^V“^s“^ 
were obtained. Metallicity parameter kpl of LuPdBi is then equal to 5. 
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FIG. 8. Analysis of the SdH oscillations in LuPdBi. (a) Lifshitz-Kosevich description (solid 
line) of the temperature dependence of the SdH amplitude. Inset: Dingle plot for the oscillations 
measured at 2.5 K. Solid straight line determines the Dingle temperature (see the text), (b) 
Landau-level fan diagram for the SdH oscillations measured at 2.5 K. Fitted straight solid line 
intersecting the ordinate axis at 0.4 indicates the Berry phase very close to tt. Dashed (blue) line 
represents the ht to extrema in conductivity AcJxx derived with macroscopic Hall coefficient, i?H- 
Dot-dashed (red) line corresponds to Auxx derived with Hall coefficient of surface channel Rs (see 
text for explanation). Inset shows the blow up of the area where the fitted lines intercept ordinate 
axis. 

Figure presents the Landau-level fan diagram for the SdH oscillations in LuPdBi at 
2.5 K. Positions of the extrema in Ap^{l/B) are marked with circles as a function of their 
number (n for maxima and n + ^/2 for minima). Though the very small value of /sdn allowed 
to observe only hve such extrema, a good linear £t could be performed, with an intercept of 
0.391(26), shown as a solid black line. However, the choice of resistance pxx or conductance 
(Txx for Landau level indexing may be very important and yield different phase factor .1^ Since 
the Hall resistivity of LuPdBi measured at 2.5 K was Pxy{B = 9T) = 0.54 mGcm, that is 
very close to Pxx = 0.50mf2cm, it cannot be ignored in analysis of the SdH oscillations. We 
thus calculated (Txx(H) = Pxx(-B)/(pL(-®) +Pxy(-®)) extracted the oscillatory component 
Aa^xi^/B). These data are shown as blue points in Figure [^. Straight line htted to the 
positions of their extrema (shown as dashed blue line in Figure]^) yielded the phase factor 
of 0.285(40) (for clarity the datapoints are not shown). 

But the case of LuPdBi is an example (which very likely can be extended to other bulk 
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samples with topologically nontrivial states on their surface) where the oscillatory component 
of o'xx(l/-B) should be derived in a different manner. Our sample, as the temperature 
dependence of resistivity indicates, contains two parallel conduction channels: metallic and 
semiconducting (Figure [^). At low temperatures, the latter makes negligible contribution 
to the conductivity, but may still play an important role in the Hall effect. However, the SdH 
oscillations we observe seem to originate from surface carriers. Hence, the Hall resistivity 
of the surface channel should be used for Acrxx(l/-B) calculation. It can be easily obtained 
as Ps^xy{B) = RgB, where Rg is the Hall coefficient of the surface channel, which can be 
derived from the formula Rg = t/{en 2 B), where n 2 D is the surface carrier concentration and 
t is the sample thickness.® For our sample of LuPdBi,/?^ = 0.0146 mf2cm/T at 2.5 K. After 
recalculating as p^^iB)/+ p^^^^^B)), extracting its oscillatory component 

(shown with red in Figure]^) and hnding positions of its extrema onl/B axis, we could £t 
another straight line resulting in the phase factor of 0.364(44) (note dot-dashed red line in 
Figure |^). 

Most remarkably, all the applied methods of determinig the phase factor (we believe that 
the electronic nature of 3D TIs is most appropriately accounted for in that employing Rg) 
yielded values very close to ^/ 2 , which was theoretically predicted for Dirac fermions. Small 
departure of the experimental value from Y 2 may arise from Zeeman coupling effedP^^ or 
small deviation of the energy dispersion of Dirac fermions from the linear one.^^ Nevertheless, 
the observation in LuPdBi of a tt Berry phase is another strong indicator of the topologically 
nontrivial character of the compound. 


CONCLUSIONS 

We synthesized high-quality single crystals and investigated the electronic properties of a 
putative 3D topological superconductor LuPdBi. The electrical transport in the compound 
was found to be driven by two parallel conducting channels: metallic and semiconducting, 
with the negligible contribution of the latter at low temperatures. The superconductivity 
emerges at R = 1.9 K and is characterized by the upper critical held i?c2(0) < 2.3 T. It is 
of type H clean-limit BCS. Remarkably, no feature in the heat capacity is observed near the 
onset of superconducting state, which rather rules out its bulk character. Instead, we tend 
to attribute the superconductivity in LuPdBi to the surface states. 
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The compound shows the linear-in-i? magnetoresistance, as high as 210% at 2.5 K and 
over 30% at 300 K (in magnetic held of 9 T). In low magnetic helds, there clearly occurs 
the weak antilocalization effect. Both features are hallmarks of the topologically protected 
states. Furthermore, the observed SdH oscillations yield the very small effective mass of 
0.06 rrie and the Berry phase very close to vr, in concert with the topological nature of 
LuPdBi, characterized by the Fermi surface sheet containing massless, extremely mobile 
Dirac fermions. 

In summary, based on the results of our study, we postulate that at low temperatures, 
below about 10 K, LuPdBi effectively becomes a topological insulator, with negligible con¬ 
tribution of its bulk to the electronic transport and with the topological nature of its surface 
states being clearly rehected in the physical properties. The superconductivity in LuPdBi 
seems to emerge from the topologicaly protected surface states and further conhrmation, 
both experimental and theoretical, of this conjecture is a challenging task for the future. 


METHODS 

Material preparation 

Single crystals of LuPdBi were grown from Bi flux. The elemental constituents of high 
purity (Lu: 99.99 wt.%, Pd: 99.999 wt.%, Bi; 99.9999 wt.%), taken in atomic ratio 1:1:12, 
were placed in an alumina crucible, and sealed under argon atmosphere inside a molybdenum 
ampule. The reactor was heated slowly up to 1300 °C and kept at this temperature for 15 
hours. Next, it was cooled down to 1000 °C at a rate of 3°C/hour, and then to room 
temperature at a rate of 5°C/hour. The crucible with the charge was extracted from the 
molybdenum ampule, covered with another alumina crucible hlled with silica wool and 
resealed in an evacuated quartz tube. Then, the tube was heated up to 900 K at which 
point the excess of Bi flux was removed by centrifugation. The crystals obtained by this 
technique had a shape of cubes or plates, with dimensions up to 0.5 x 0.5 x 0.5 mm^ or 
0.1 X 0.5 X 3mm^, respectively. They had metallic luster and were stable against air and 
moisture. 
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Material characterization 


Powdered single crystals of LuPdBi were characterized at room temperature by powder 
X-ray diffraction (PXRD) carried out using an X’pert Pro PANanalitical diffractometer with 
Cu-Ka radiation. The crystal structure refinements and the theoretical PXRD pattern cal¬ 
culations were done employing the FULLPROF program.^ Quality of the samples studied 
was also verified by single crystal X-ray diffraction (SCXRD) performed at room temper¬ 
ature on an Oxford Diffraction Xcalibur four-circle diffractometer equipped with a CCD 
camera and using Mo-Ka radiation. The SCXRD data analysis was done using the program 
package SHELXL-97.® The results of PXRD and SCXRD investigations are presented in 
Supplementary Material. 

Chemical composition of the LuPdBi crystals was examined on a FEI scanning electron 
microscope (SEM) equipped with an EDAX Genesis XM4 spectrometer. The specimens 
were glued to a SEM stub using carbon tape. The crystals were found homogeneous and 
free of foreign phases, with the chemical compositions very close to the ideal one (for details 
see Supplementary Material). 


Physical measurements 

Electrical resistivity and Hall measurements were carried out from 0.4 K to 300 K in 
applied magnetic fields up to 9 T using a conventional ac four-point technique implemented in 
a Quantum Design PPMS platform. Dimensions of samples were: 0.010 x 0.051 x 0.227 cm^, 
for pxx, and 0.0067 x 0.051 x 0.1524 cm^ for pxy measurements. Current and voltage leads 
were 50 pm thick silver wires attached to the parallelepiped-shaped specimens with silver 
paste and additionally spot welded. The heat capacity was measured on 10.6 mg collection 
of single crystals, in the temperature interval 0.4—5 K by relaxation time method using 
also the PPMS platform. DC magnetization and AC magnetic susceptibility measurements 
were performed on 30.1 mg collection of single crystals, in the temperature range 1.62—2.6 
K in weak magnetic fields up to 2 mT employing a Quantum Design MPMS-XL SQUID 
magnetometer. 
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SUPPLEMENTARY MATERIAL 


MATERIAL CHARACTERIZATION 


X-ray diffraction 


The PXRD results (see Fig. SI) confirmed a single-phase character of the obtained single 
crystals of LuPdBi. The X-ray diffraction pattern can be fully indexed within the E43m 
space group, characteristic of half-Heusler compounds, and yields the cubic lattice parameter 
a = 6.565(1) A. This value is in perfect accord with the literature value 6.566(1) A deter¬ 
mined for poly crystalline sample,^ and the value 6.56 A reported for thin dims of LuPdBi.® 
It is however markedly different from the lattice parameter of 6.63 A stated by Xu et ah for 
their powdered single crystals.® As can be inferred from Fig. SI, the experimental PXRD 
pattern of our single-crystalline LuPdBi can be very well modeled with the MgAgAs-type 
crystal structure with the Lu atoms located at the crystallographic 4a (0,0,0) sites, the Pd 
atoms occupying the 4c (1/4,1/4,1/4) sites, and the Bi atoms placed at the 46 (1/2,1/2,1/2) 
sites. The same structural model was recently shown to account well for the PXRD data of 
powdered single crystals of YPdBi and YPtBi.® 



Figure SI. XRD pattern for powdered single crystals of LnPdBi. 
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Table SI. Crystal data and structure refinement for LuPdBi. 


Fonnula weight 
Temperature 
Wavelength 
Crystal system 
Space group 
Unit cell parameter 
Volume 
Z 

Density (calculated) 

Absorption coefficient 
F(OOO) 

Crystal size 

Theta range for data collection 
Index ranges 
Reflections collected 
Independent reflections 
Completeness to theta = 25.242° 
Absorption coiTection 
Max. and min. transmission 
Refinement method 
Data restraints parameters 
Goodness-of-fit on F- 
Final R indices [I>2sigma(I)] 

R indices (all data) 

Absolute stiucture parameter 
Extinction coefficient 
Largest diff peak and hole 


490.35 
295(2) K 
0.71073 A 
Cubic 
F43m 

a = 6.5683(10) A 
283.37(13) A^ 

4 

11.494 Mg in^ 

102.423 imn-i 
800 

0.04 X 0.03 X 0.03 imn^ 

5.377 to 26.497°. 

-8<=/7<=8, -8<=/c<=8, -8<=/<=8 
627 

48 [R(int) = 0.0808] 

96.0 % 

Numerical 

0.137 and 0.060 

Full-matrix least-squares on F- 

48 / 0'5 

1.047 

Ri = 1.76%, wR, = 4.53% 

Ri = 1.76%, wR2 = 4.53% 
-0.043(16) 

0.0048(7) 

1.079 and -1.524 e/A^ 


Table S2. Atomic coordinates and equivalent isotropic displacement parameters U{eq) for 
LuPdBi. U{eq) is defined as one third of the trace of the orthogonalized f/*-^ tensor. 


Atom 

.Y 



L (eq) [A^xlO^] 

Lu 

0 

0 

0 

9(2) 

Pd 

0.25 

0.25 

0.25 

10(2) 

Bi 

0.5 

0.5 

0.5 

7(1) 
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It is worth noting that despite the aforementioned divergence in the values of the lattice 
parameter a, the X-ray diffractogram of LuPdBi presented in Ref. El appears very similar to 
our result displayed in Fig. SI. Small differences in the relative intensities of a few Bragg 
peaks [e.g., (311), (222), (420)] might be attributed to a certain level of atomic disorder in 
the sample studied by Xu at ah, since our PXRD data were rehned to a small value of the 
residual Bragg factor R = 4.61%, with full occupancies of all the crystallographic positions 
and no structural disorder (see also the discussion presented in Ref. Ej). 

The crystal structure of the studied samples of LuPdBi was also checked by the single 
crystal X-ray diffraction method. Details on the performed rehnement are gathered in Table 
SI. The positional data and the equivalent isotropic displacement parameters are given in 
Table S2. The so-re£ned lattice parameter a = 6.5683(10) Ais very close to that derived from 
the PXRD data, and thus signihcantly smaller than the result reported for single-crystalline 
LuPdBi by Xu at al.^ 


Energy-dispersive X-ray spectroscopy 

Examples of the EDX spectrum and the SEM image obtained for the single crystals 
of LuPdBi is presented in Fig. S2. The chemical composition derived from these data is 
Lu 32 ,oo( 32 )Pd 33 . 93 ( 33 )Bi 34 07 ( 34 ), in a Very good accord with the ideal equiatomic one. 



Figure S2. SEM-EDX results for single-crystalline LuPdBi. 
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ELECTRONIC AND THERMODYNAMIC PROPERTIES 


WAL parameters 

2 

The temperature dependencies of the weak antilocalization parameters and (to- 

2 

gether with their uncertainties d{^^) and dL^p) obtained from fitting the Hikami-Larkin- 

Nagaoka (HLN) formula (Eq. 1 in the main text) to the experimental u{B) data of LuPdBi 

are collected in Table S3. When both parameters were free to vary in htting procedure, 

the values given in second and fourth columns of the Table S3 were obtained. However the 

varied with temperature in irregular way and both parameters were strongly dependent 

on each other, as shown by statistical dependency values in fifth column. Thus we could 

2 

not consider all these values of and as reliable. Reasoning that the WAL model is 

most adequate at lowest temperature, where the contribution from the bulk semiconducting 

2 

channel is the smallest, we repeated the htting with the value of hxed at the value 
derived at T = 2.5 K. Resulting values of are collected in sixth column of the Table. It 
is worth noting that their errors are now smaller by about one order of magnitude. These 
values decrease monotonously with increasing temperature, which seems to rehect the in¬ 
creasing of the bulk channel contribution. 

TABLE S3. WAL parameters for sheet conductance of LuPdBi at different tempera¬ 
tures. 


T(K) 

rie^ 

2'K‘^h 


Lip 

dLp 

dependency 

Lip\ri{T=2.5K) 

dL^\jj(^T=2.5K) 

2.5 

-153.7 

21.8 

49.6 

2.2 

0.9976 

49.6 

0.1 

4 

-340.5 

149.1 

33.1 

4.0 

0.9996 

41.3 

0.1 

7 

-99.8 

21.6 

43.5 

2.8 

0.9986 

38.4 

0.1 

10 

-65.6 

12.9 

48.8 

3.0 

0.9977 

37.9 

0.1 

50 

-49.6 

7.4 

44.7 

2.0 

0.9984 

32.6 

0.1 

150 

1 

cc 

bo 

0.5 

52.1 

0.9 

0.9971 

23.6 

0.1 
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Figure S3. Temperature variations of the WAL parameters for single-crystalline LuPdBi. Green 

2 

triangles represent values of circles and diamonds correspond to obtained from two 

fitting procedures described in text. 


2 

For more clarity we plotted the values of and L^f, resulting from both htting proce¬ 
dures, versus temperature, as shown in Fig. S3. 

The temperature variation of the phase coherence length Lip\^(T= 2 .bK) (resulting from the 
2 

HLN £t with fixed at the value derived at T = 2.5 K.) can be roughly described by the 
formula given in Ref. [3] (represented in Fig. S3 by the solid red curve), however such a fit 
yields the electron-phonon interaction parameter close to zero. 

The parameter derived in Ref. |3|was about an order of magnitude smaller than that 
obtained on our sample, whereas corresponding values of were about twice as large as our 
L^\ri{T= 2 . 5 K) data. This is rather a good agreement, taking into account different electronic 
character of the samples examined in our work and that described in Ref. El It seems to be 
a common feature of half-Heusler compounds that the values of parameter —r] are of order 
10® — 10®, thus much larger than 1/2 expected for 2D TI systems. Apart from LuPdBi it 
has very recently been observed in samples of LuPtSb,® YPtBi and LuPtBi.® Authors of all 
these papers attributed this fairly unexpected result to contributions from bulk or side wall 
conductivity channels. 
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Hall parameters 


In concert with the previous report,^ the Hall carrier concentration in the single crys¬ 
tal investigated in our work was found temperature independent below 50 K and then to 
increase with increasing temperatures. However, in the entire temperature range covered, 
the magnitude of tt-h determined in our study is approximately four times smaller than that 
estimated by Xu et ah (their result was nn — 4.8 x 10^®cm“^ at low temperatures). Re¬ 
markably, our hndings are completely different as regards the Hall carrier mobility. For our 
sample, /th was found to gradually decrease with increasing temperature from the value 
2404cm^V“^s“^ at T = 2.5 K to 573cm^V“^s“^ at 300 K, while Xu et ah reported an 
increase from 330 to 380 cm^V“^s“^ in the same temperature interval. 


Heat capacity 

The striking result reported in Ref. |3]is the occurrence of a pronounced A-shaped anomaly 
in the specihc heat of LuPdBi that coincides with the superconducting transition. Conse¬ 
quently, Xu et ah ascribed the superconductivity in their single crystals to the bulk. This 
result markedly differs from the behavior of our single crystals, where no visible feature in 
C{T) was found near Tc, despite the formation of superconducting state was convincingly 
proved by means of the magnetic and electrical transport measurements. This finding led 
us to the conclusion that the superconductivity in our samples of LuPdBi is confined to the 
surface, with the superconducting condensate volume negligibly small as compared to the 
bulk. The result obtained by Xu et ah differs also from the properties of closely related 
systems LuPtBi and YPtBi,® for which the specific heat was found featureless at the su¬ 
perconducting transition. For the latter compound we recently conhrmed the lack of any 
anomaly at Tc in our own investigation carried out on high-quality single crystals.® 

The observed discrepancy might be related to the different electronic character of the 
single crystals studied in our work and those investigated in Ref. El Comparison of the 
temperature-dependent electrical resistivity suggests that our samples are much more 
semimetallic than the predominantly semiconducting one reported by Xu et ah In this 
context, however, the observed over-an-order-of-magnitude discrepancy between the Som- 
merfeld coefficient 7 = 0.75mJ/molK^ derived in the present work and 7 = 11.9mJ/molK^ 


found by Xu et al. becomes most surprising. Here, again, one should note that the electronic 
specihc heat of the related semimetallic nonmagnetic bismuthides, like LuPtBi and YPtBi 
(Refs. I7 HTU]) . and antimonides, like LuPtSb, YPdSb and YPtSb (Ref. |TT]) is similar to that 
established by us for LuPdBi. 
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